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We report a neutron scattering study of the spin correlations for the spin 5/2, two-dimensional 
antiferromagnet Rb2MnF4 in an external magnetic field. Choosing fields near the system's bicritical 
point, we tune the effective anisotropy in the spin interaction to zero, constructing an ideal S = 5/2 
Heisenberg system. The correlation length and structure factor amplitude are closely described by 
the semiclassical theory of Cuccoli et al. over a broad temperature range but show no indication 
of approaching the low-temperature renormalized classical regime of the quantum non-linear sigma 
model. 
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Magnetic systems with reduced dimensionality have 
provided a basis for numerous insights into the varying 
roles quantum and thermal fluctuations play in driving 
phase transitions. The two dimensional quantum Heisen- 
berg antiferromagnet (2DQHA) represents a particularly 
important example of such a system due to the possible 
connections between it and the origins of superconduc- 
tivity in the layered copper oxides 0. Specifically, the 
parent compounds to the high T c superconductors, such 
as La2Cu04, are good examples of 2DQHA's on a square 
lattice with spin, S = 1/2, and the magnetism persists 
into the superconducting state |Q. As a consequence of 
this connection, focus on the 2DQHA has intensified in 
recent years, and the combined efforts of experiment, the- 
ory, and simulation have yielded an increasingly cohesive 
picture of its behavior. 

In mapping the 2DQHA onto the quantum nonlinear 
sigma model (QNLcM), Chakravarty, Halperin, and Nel- 
son have produced an effective field theory for the sys- 
tem which predicts renormalized classical behavior at 
low temperatures. Calculations from this theory (la- 
belled CHN-HN Q), describe a quasiexponential growth 
in the correlation length, £, with inverse temperature 
that matches measurements on S — 1/2 systems Q with 
good absolute agreement over a remarkably large tem- 
perature range. However, high temperature results from 
systems with S > 1/2 vary markedly from the CHN- 
HN predictions ^-|). These discrepancies have led to 
questions regarding the applicability of the effective field 
theory in the temperature range probed in these experi- 
ments and have motivated new approaches to describing 
the 2DQHA at high temperatures. For example, Elst- 
ner et al. JtJ have applied series expansion techniques 
to characterize the onset of correlations with good accu- 
racy. In addition, Cuccoli et al. [|| have recently proposed 
a semiclassical theory, the pure quantum self consistent 
harmonic approximation (PQSCHA). This theory calcu- 



lates thermodynamic quantities, such as £, with respect 
to their classical values, with the effects of quantum fluc- 
tuations entering as a renormalization of temperature. 
While this semiclassical picture should be valid only at 
high temperature for S — 1/2, the authors predict that 
the theory should work over a broad range of tempera- 
tures for S > 1/2 §. 

A major obstacle to forming a unified picture of the 
2DQHA, in which the relative merits of these theoretical 
approaches can be distinguished, has been the absence 
of experimental systems with S > 1/2 that model the 
2DQHA over a broad range of length scales and, con- 
comitantly, temperatures. In particular, the materials 
that have been studied, which include K2MF4 || and 
La 2 Ni0 4 (§ with S = 1 and Rb 2 MnF 4 §| and KFeF 4 § 
with S = 5/2, possess appreciable Ising anisotropies that 
generate crossovers from 2D Heisenberg to 2D Ising be- 
havior at quite short length scales. For example, in both 
Rb2MnF 4 and KFeF 4 the crossover to Ising behavior is 
first apparent at £/a ~ 6, where a is the lattice spac- 
ing. In contrast, the 2D S — 1/2 material Sr2Cu02Cl2 
exhibits model 2DQHA behavior in its spin correlations 
for length scales up to £/a ~ 200 [||. As spin number is 
the central parameter controlling the strength of quan- 
tum effects, measurements over a broad range of £ away 
from the extreme quantum limit, S = 1/2, are crucial to 
enhance our understanding of the 2DQHA. 

In this paper we present a study of Rb2MnF 4 under 
conditions in which its spin space anisotropy is effec- 
tively reduced to zero, thus allowing us to characterize 
the Heisenberg behavior of this 5 = 5/2 system over a 
range of length scales which rivals those probed in the 
experiments on the S — 1/2 systems. Rb2MnF 4 has the 
K2NiF 4 crystal structure, in which magnetic planes of 
MnF2 form a square lattice with nearest neighbor spac- 
ing, a = 4.2152 at 5 K. The Mn 2 + ions (S = 5/2) inter- 
act antiferromagnetically with nearest neighbors on the 
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lattice with an isotropic exchange, J = 7.36 K, plus a 
uniaxial, dipolar field along the c axis (perpendicular to 
the planes) of strength, q^bHa ~ 0.005 J |J. 

To reduce the effective anisotropy we place Rb2MnF4 
in an external magnetic field, H, parallel to c. As stud- 
ied in detail by Cowley et al. |l(J, the phase diagram for 
Rb2MnF4 in a uniform magnetic field contains a spin- 
flop transition between the low-field Ising phase and a 
high-field XY phase. Results from Ref. jl0| along with 
new measurements, shown in Figure 1, reveal the T = 0, 
2D bicritical point, where the boundaries between these 
phases and the paramagnetic phase meet, at Hb = 
5.30 ± 0.12 T. The effective anisotropy, g, in the spin 
interaction varies with field as g ~ H 2 — H B — CT jio) , 
with g < denoting easy axis anisotropy and g > de- 
noting easy plane. Thus, tuning the field at a given tem- 
perature provides a powerful mechanism for controlling 
the symmetry and magnitude of the effective anisotropy. 
In particular, for the fields at which g — (the dashed 
line in Fig. 1), the effective anisotropy vanishes, and all 
spin components are critical. In our experiment we have 
carefully chosen fields to match this condition as closely 
as possible, thus avoiding a crossover to a lower symme- 
try transition. Specifically, we obtain H B and C from 
the spin-flop line at low temperature (T < 25 K) and 
extrapolate this line into the paramagnetic region to de- 
termine the appropriate fields. This strategy has allowed 
us to track the Heisenberg behavior in Rb2MnF4 to sig- 
nificantly lower temperatures and, concomitantly, larger 
length scales than otherwise possible. 

We have characterized the instantaneous magnetic 
correlations of this Heisenberg system by performing 
2-axis, energy-integrating neutron scattering measure- 
ments with the BT9 spectrometer at the NIST Center 
for Neutron Research. To accomodate the relative orien- 
tations of the spectrometer axes and the external field, 
we oriented the crystal such that c was perpendicular to 
the scattering plane. To guarantee that the quasielastic 
approximation was satisfied, we employed a large inci- 
dent neutron energy, Ei n — 100 meV. We created the 
100 meV incident beam using the (0,0,4) reflection from 
a pyrolytic graphite monochromator. The collimator se- 
quence for the spectrometer was 40'-13'-S-10', and a sap- 
phire filter after the monochromator removed higher or- 
der contamination. We model the static structure factor, 
S i (Q2d)) where Q2D is the momentum transfer in the 
magnetic planes, with a simple Lorentzian lineshape 

J(Q 2 d) oc S(Q 2D ) = ^° , (1) 

where q 2 D is the wave vector measured from the antifer- 
romagnetic zone center, q2D = Q2D — (1/2, 1/2,0). Fig- 
ure 2 displays scattering profiles from the 2-axis scans 
along the direction (K/2, — K/2, 0) at three representa- 
tive temperatures. We fit the scattering intensity to Eq. 



1 convolved with the spectrometer resolution function to- 
gether with a sloping background. The fit results, shown 
with solid lines in Fig. 2, describe the data very well 
at all temperatures. In addition, to test the validity of 
the quasielastic approximation, we have characterized the 
shape of the dynamic structure factor, S(Q2D,w), in en- 
ergy transfer, u>, by performing additional 3-axis mea- 
surements scanning lo at q 2 D = 0. A Lorentzian line- 
shape in uj fits these 3-axis scans well. An analysis of 
the 2-axis scans which accounts explicitly for the shape 
of the dynamic structure factor gives results for £ and So 
that agree with those obtained assuming the quasielas- 
tic approximation to within 5% at all temperatures, thus 
confirming the validity of the approximation for our mea- 
surements. 

Figure 3a shows £/a extracted from the fits on a log 
scale versus JS(S + 1)/T. We choose JS(S + 1) rather 
than JS 2 for the energy scale following the observation 
by Elstner et al. [0 that £ depends on S primarily through 
JS(S + 1) at high temperature for S > 3/2. Also plotted 
are the results from Lee et al. 0] for Rb2MnF4 in zero 
magnetic field. The good agreement between the two 
measurements in the temperature region in which they 
overlap confirms that our method to determine the in- 
stantaneous spin correlations is accurate. In addition, as 
the plot illustrates, by reducing the effective anisotropy 
to zero we have been able to extend significantly the tem- 
perature range of the Heisenberg behavior. At our lowest 
temperature £/a exceeds 100, rivaling the experimental 
range for the best 5 = 1/2 Heisenberg systems. 

The solid line in Fig. 3a is the prediction for £/a from 
the PQSCHA ||. As the figure demonstrates, the cal- 
culated behavior from the semiclassical theory agrees 
closely with the measured correlation lengths over the full 
experimental range. As mentioned above, the PQSCHA 
calculates £ with respect to its value for the classical 
Heisenberg model. For reference, the dashed-dotted line 
shows the correlation length for the classical model, £c l > 
as obtained from simulation [Q|ll]]. At high temperature 
(£/a < 5), £ closely approximates £cl when the energy 
scale is taken as JS(S+ 1). Sokol et al. |L2l have labelled 
this temperature region the 'classical scaling regime'. 
However, as the figure demonstrates, with decreasing 
temperature quantum effects become increasingly impor- 
tant, and this simple S(S + 1) scaling breaks down. The 
renormalization of temperature that the PQSCHA pro- 
vides accounts for the quantum effects nicely. 

Sokol et al. have argued that the 5 = 5/2 2DQHA 
should cross over from classical scaling to renormalized 
classical behavior near T ~ A, where A is the upper 
bound of the spin wave energy spectrum (approximately 
5.5 meV at 35 K). Thus, we estimate this temperature 
for our system to be near 60 K (JS(S + 1)/T « 1.1), 
roughly where the deviations between £cl and the ex- 
perimental results in Fig. 3a become noticeable. How- 
ever, the low-temperature, renormalized classical predic- 
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tion for the QNLcrM from CHN-HN g, shown with the 
dashed line in Fig. 3a, deviates significantly from the 
measured correlation lengths even for ^/a exceeding 10 2 . 
This result contrasts with those for the 5 = 1/2 systems 
which show good agreement with the theory over length 
scales down to £/a ~ 1. However, recent quantum Monte 
Carlo simulations combined with finite size scaling have 
indicated that this experimental agreement is in part ac- 
cidental due to a cancellation between finite lattice ef- 
fects and higher order corrections to the QNLcrM Jl3| . 
In addition, Elstner et al. have argued that the onset 
of the renormalized classsical behavior described by the 
QNLcrM should occur at progressively higher £ for larger 
5 JtJ. Our measurements on this 5 = 5/2 2DQHA, show- 
ing deviations from the projected asymptotic behavior 
even as the measurements extend well below the classical 
scaling regime, support the view of an extended crossover 
region between the classical and renormalized classical 
scaling regimes. 

In Figure 3b we display Sq, the amplitude of the static 
structure factor at the ordering wave vector, extracted 
from the measurements. The zero field results [El, shown 
with open triangles, have been placed on an absolute 
scale by setting the high temperature limit to 5(5+ 1)/3. 
Our results, placed on the same scale, extend the data 
smoothly to lower temperature. The solid line in Fig. 3b 
is the prediction from the PQSCHA with no adjustable 
parameters. As with £, the PQSCHA predicts So accu- 
rately over a broad temperature range. 

For both the classical Q and quantum 2DQHA's, 
low-temperature theories predict 5o/£ 2 oc T 2 . How- 
ever, experiments over the range 0.2 < T/2np s < 0.4 
for 5 = 1/2 and 5 = 1 (§ have found 5 /£ 2 to de- 
pend weakly on T. The inset to Fig. 3b shows the results 
for this ratio for Rb2MnF4. As the figure illustrates, at 
high temperature 5o/£ 2 is roughly constant, but at lower 
temperature the ratio has a strong temperature depen- 
dence, qualitatively consistent with the low-temperature 
theories. The break from temperature independent 5o/£ 2 
occurs near T/2Tip s — 0.2, the lower limit of results on 
other experimental systems. Thus, the measurements to 
low temperature (T/2irp s ~ 0.11) in Rb 2 MnF 4 appear 
to have captured a crossover from this high temperature 
regime. 

Finally, we note that the bicritical phase boundaries 
themselves should reflect the 2D correlations. Specifi- 
cally, the Ising and XY critical lines should approach 
each other tangentially with g ~ T~ 2 £~ 2 JTo| . The solid 
lines in Fig. 1 correspond to this form with £ from Fig. 
3a. In the fit the g = line is fixed by the low temper- 
ature spin-flop data, and only the Ising and XY ampli- 
tudes are adjusted. Clearly, the agreement is good for 
temperatures below ~ 35K. 

In conclusion, the presence of a bicritical point in 
the phase diagram of Rb2MnF4 has provided a direct 
and convenient mechanism for controlling the effective 



anisotropy in the spin interaction of this two-dimensional 
5 = 5/2 antiferromagnet. In the present study we have 
focused on conditions in which the effective anisotropy is 
essentially zero, allowing us to track the growth of cor- 
relations in this large-spin, two-dimensional Heisenberg 
system to correlation lengths exceeding 100 lattice units. 
The success of this approach, the first to probe such ex- 
tended Heisenberg correlations in a system with 5 > 1/2, 
suggests that similar studies near the bicritical point in 
tetragonal systems with other spin quantum numbers 
would be fruitful in further extending our understanding 
of the 2DQHA. Our measurements for 5 = 5/2 demon- 
strate the success of a semiclassical theory, the PQSCHA, 
over a broad temperature range. The disagreement be- 
tween the measured behavior and the predictions of the 
effective field theory based on the QNLcrM, even as £ be- 
comes large, confirms speculations that the temperature 
range applicable to this theory lies well below that ac- 
cessible in experiment. However, the growing departure 
from the simple scaling to classical behavior in Fig. 3a 
along with the marked temperature dependence of 5o/£ 2 
suggests behavior distinct from that seen at high temper- 
ature. Thus, the experiment reveals an extensive temper- 
ature region for the crossover between the classical and 
renormalized classical regimes. A quantitative theory for 
this crossover would be most valuable. 
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FIG. 1. Phase diagram for RD2M11F4 in a magnetic field, 
H, parallel to the c axis. The solid circles (shifted by 0.15 
T) are from Cowley et al. [10], and the open circles are from 
the current experiment. The solid lines are a fit to the phase 
boundaries, g ~ T~ £~ , using the results for £ from Fig. 
3a. Measurements of the spin correlations at zero effective 
anisotropy are taken along the dashed line, which extends 
from the T = 0, 2D bicritical point at H B = 5.30 ± 0.12 T. 
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FIG. 2. 2-axis scans along the direction (K/2,-K/2,0) at 
field and temperature values for which the effective anisotropy 
in the spin interaction is near zero. The solid lines are fits to 
Eq. 1 convolved with the instrumental resolution function. 
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FIG. 3. (a) The magnetic correlation length, £, in units 
of the lattice spacing, a, from the present experiment (solid 
circles) and measured at zero field in the traditional scat- 
tering geometry (open triangles) [5]. The solid line is the 
prediction from the semiclassical theory of Cuccoli et al. [8]. 
The dashed-dotted line is the correlation length for the clas- 
sical Heisenberg model obtained from simulation [7,11]. The 
dashed line is the prediction from low-temperature QNLctM 
effective field theory of CHN-HN [2]. (b) So, the ampli- 
tude of the static structure factor at the ordering wave vector, 
from the present experiment (solid circles) and measured at 
zero field in the traditional scattering geometry (open trian- 
gles) [5] . The solid line is the prediction from the semiclassical 
theory of Cuccoli et al. [8]. The inset shows So/£ 2 versus tem- 
perature. The dashed line is a fit to the solid circles assuming 

s /e ~ t 2 . 
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